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ABSTRACT 



We study the amphfication of electromagnetic vacuum fluctuations induced by the evolution 
of scalar metric perturbations at the end of inflation. Such perturbations break the conformal 
invariance of Maxwell equations in Friedmann-Robertson- Walker backgrounds and allow the 
growth of magnetic fields on super-Hubble scales. We relate the strength of the fields generated 
by this mechanism with the power spectrum of scalar perturbations and estimate the amplifica- 
tion on galactic scales for different values of the spectral index. Finally we discuss the possible 
effects of finite conductivity during reheating. 



1 Introduction 



The existence of cosmic magnetic fields with large coherence lengths (> 10 kpc) and typical 
strength of 10~^ G, still remains an open problem in astrophysics [|I|. A partial explanation, 
widely considered in the literature, is based on the amplification of seed fields by means of 
the so called galactic dynamo mechanism. In this mechanism, the differential rotation of the 
galaxy is able to transfer energy into the magnetic field, but nevertheless it still requires a 
pre-exisiting field to be amplified. The present bounds on the necessary seed fields to comply 
with observations are in the range Bgeed ^ 10~^^ — 10~^^ G {h = 0.65 — 0.5) at decoupling 
time, coherent on a comoving scale of ~ 10 kpc, for a flat universe without cosmological 
constant. For a flat universe with nonvanishing cosmological constant, the limits can be relaxed 
up to Bseed > 10^2^ - 10-30 G {h = 0.65 - 0.5) at decoupling for = 0.7 and fi^ = 0.3 §. 
The observations of micro-Gauss magnetic fields in two high-redshift objects (see [I], |^ and 
references therein) could, if correct, impose more stringent conditions on the seeds fields or 
even on the dynamo mechanism itself. 

The cosmological origin of the seed fields is one of the most interesting possibilities, although 
some other mechanisms at the astrophysical level, such as the Biermann battery process, have 
also been considered 0, In the cosmological case, in which we will be mainly interested 
in this work, it is natural to expect |^ that the same mechanism that gave rise to the large- 
scale galactic structure, i.e. amplification of quantum fluctuations during inflation, was also 
responsible for the generation of the primordial magnetic fields. However, it was soon noticed 
that the gravitational amplification does not operate in the case of electromagnetic (EM) 
fields. This is because of the conformal triviality of Maxwell equations in Friedmann- Robertson- 
Walker (FRW) backgrounds, i.e. conformally invariant equations in a conformally flat space- 
time. In order to avoid this difficulty, several production mechanisms have been proposed in 
which Maxwell equations are modified in different ways. Thus for example, the addition of 
mass terms to the photon or higher-curvature terms in the Lagrangian was studied in 0. The 
contribution of the conformal anomaly was included in . In the context of string cosmology, 
the effects of a dynamical dilaton field were taken into account in 0. Other examples include 
non-minimal gravitational-electromagnetic coupling [^], inflaton coupling to EM Lagrangian 0, 
spontaneous breaking of Lorentz invariance [|10| or backreaction of minimally coupled charged 



scalars ||TT], |T2|, 0. Some of them are able to generate fields of the required strength to seed 



the galactic dynamo or even to account for the observations without further amplification. 

In this paper we explore the alternative possibility, i.e. we avoid conformal triviality by con- 
sidering deviations from the FRW metric (see for a suggestion along these lines). This ap- 



proach is rather natural since we know that galaxies formed from small metric inhomogeneities 
present at large scales and, in addition, it does not require any modification of Maxwell electro- 
magnetism. In the inflationary cosmology, metric perturbations are generated when quantum 
fluctuations become super-Hubble sized and thereafter evolve as classical fluctuations, reenter- 
ing the horizon during radiation or matter dominated eras The same mechanism would 



operate on large-scale EM fluctuations. However, if conformal invariance is not broken, each 
positive or negative frequency EM mode will evolve independently, without mixing. This im- 
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plies that photons cannot be created and therefore magnetic fields are not amplified. However, 
in the presence of an inhomogeneous background, we will show that the mode-mode coupling 
between EM and metric perturbations generates the mixing. This in turn will allow us to relate 
the strength of the magnetic field created by this mechanism and the particular form of the 
metric perturbations described by the corresponding power spectrum. Those photons produced 
in the infiation-radiation transition with very long wavelengths can be seen as static electric 
or magnetic fields. Because of the high conductivity of the Universe in the radiation era, the 
electric components are rapidly damped whereas, thanks to magnetic fiux conservation, the 
magnetic fields will remain frozen in the plasma and their subsequent evolution will be trivial, 
Ba? =const P, |^ . The paper is organized as follows. In section 2 we obtain the Maxwell equa- 
tions in the presence of an inhomogeneous background and calculate the occupation number 
of the photons produced. In section 3 we apply these results to calculate the corresponding 
magnetic field generated at galactic scales. Section 4 is devoted to the analysis of the effects 
of finite conductivity in those results and finally, section 5 includes the main conclusions of the 
paper. 



2 Maxwell equations and photon production 



Although there are previous works on the production of scalar and fermionic particles in inho- 
mogeneous backgrounds |jl6|, [l^, in this paper we will need to extend the analysis to the case 
of gauge fields. Let us then consider Maxwell equations 

V,F>^^ = 0, (1) 

in a background metric that can be splitted as g^i, = g'^^ + h^^, where 

g^^^dx^dx" = a^{r]){d'rf - dijdx'dx^) (2) 

is the fiat FRW metric in conformal time and 

h^udx^dx" = 2a^{r])^{drf' + 5ijdx^dx^) (3) 

is the most general form of the linearized scalar metric perturbation in the longitudinal gauge 
and where it has been assumed that the spatial part of the energy-momentum tensor is diagonal, 
as indeed happens in the infiationary or perfect fiuid cosmologies |jl5|. In this expression ^{rj, x) 
is the gauge invariant gravitational potential. The equation (|I]) can be written as: 

^ ^ (^g''''g'^{d^Af,-df,A^))=0, (4) 



yfg dx^^ 

which leads in this background to the following linearized equations 



■^{{l-2<!>){d,Ao-doA,)) = 0, (5) 
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for z/ = and 



|^((l-2$)(9,Ao-9oA,)) 
+ ^ ((1 + 2$)(a,A, - d,A,)) = 0, (6) 

for V = i. In addition, we will use the Coulomb gauge condition V ■ A = 0. 

In order to study the amplification of vacuum fluctuations, let us consider a particular 
solution of the above equations that we will denote by A^^^{x) such that asymptotically in the 
past it behaves as a positive frequency plane wave with momentum k and polarization A, i.e, 

4\x) "^-^"^ Af^^\x) = ^j^e,{k, X)e^^^-'^^ (7) 

where k"^ = k"^. For the two physical polarization states we have, eik, X)-k = and eo(A;, A) = 0. 
We will work in a finite box with comoving volume V and we will take the continuum limit at 
the end of the calculation. We are assuming that metric perturbations vanish before inflation 
starts, so that we can define an appropriate initial conformal vacuum state. Because of the 
presence of the inhomogeneous background, in the asymptotic future, this solution will behave 
as a linear superposition of positive and negative frequency modes with different momenta and 
different polarizations, i.e.. 



■gv) 



It is possible to obtain an expression for the Bogolyubov coefficients akqw and Pkqw to first 
order in the metric perturbations. With that purpose, we look for solutions of the equations of 
motion in the form: 

A J'^ (x) = 4°)^'^ (x) + (x) + . . . (9) 

where A^^^'^'^^x) is the solution in the absence of perturbations given by (|^). Introducing this 
expansion in and Fourier transforming, we obtain for the temporal component of the EM 
field to first order in the perturbations: 

A^'''\q, V) = + g, ,)e--^ (10) 

where, as usual, $(g, rj) = (27r)~^/^ / (i^xe*'^^$(x, r/). The zeroth order equation implies A'Q^''''^{q, rj) 
0. The spatial equations (^ can be written to first order as: 

2<^>'Af' + 9,4')' - A^^" + 2V<^ ■ V4°) - 2V<1> ■ + V^P + 4<!>V^Af'^ = (11) 

Inserting again expansion (^, these equations can be rewritten in Fourier space as: 

^ A^'^'''(g, V) + q'A^^''\q, - jf'iq, r^) = (12) 



where: 



k"^ — k ■ q 



<l>{k + q,r]) ei{k,X)e 



~ikr) 



+ (e(^, A) ■ q) ^k + q, r^)^e-''^ - r'Stl^A. + g-, ^)e-^^) j (13) 



Solving these equations we find, up to first order in the perturbations: 

y2kV q Jm 

where rjo denotes the starting time of infiation. Comparing this expression with (|^), it is 
straightforward to obtain the Bogolyubov coefficients /Skqw, they are given by: 

P,^,y = ^== / e (g, A') ■J''\q,r^)e-^^^dr^ (15) 

where rji denotes the present time. The total number of photons created with comoving 
wavenumber kc = 2tt/Xg, corresponding to the relevant coherence length, is therefore given by 



IS 



Nka = j:T.\(3kkaXy\' (16) 
A,A' k 

We will concentrate only in the effect of super-Hubble scalar perturbations whose evolution is 
relatively simple |ll5| : 

$(fc,r/) = C,-^ (- I a'dr^+D,-^, (17) 
adr] \a J J a-^ 

the second term decreases during inflation and can soon be neglected. Thus, it will be useful to 
rewrite the perturbation as: $(/c, 77) = CkT{ri). During inflation or preheating, these perturba- 
tions evolve in time, whereas they are practically constant during radiation or matter eras. We 
will neglect the effects of the perturbations once they reenter the horizon. This is a good ap- 
proximation for modes reentering right after the end of inflation since they are rapidly damped. 
In addition, we will show that those modes are the more relevant ones in the calculation. 

The power spectrum corresponding to ( [T7| ) is given by: 

where for simplicity we have taken a power-law behaviour with spectral index n and we have set 
the normalization at the COBE scale Ac — 3000 Mpc with As — 5 ■ 10~^. In the case of a blue 
spectrum, with positive tilt (n > 1), perturbations will grow at small scales and it is necessary 
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to introduce a cut-off kmax in order to avoid excessive primordial blaclc hole production |19|. 
Accordingly, only below the cut-off the perturbative method will be reliable. For negative tilt 
or scale-invariant spectrum there will be also a small scale cut-off related to the minimum size 
of the horizon kmax ^ cliHj, where the I subscript denotes the end of inflation. 

We can obtain an explicit expression for the total number of photons (|TBp in terms of the 
power spectrum. Taking the continuum limit J2k ^ {2n)^^^'^V J d^k, we get: 



kc 



A,A' 

A,A' 



Sk 

(27r)3/2 



fcfccAA' I 



d^k \C\ 



k+kG\\ 



(27r)3/2 2kGV'^ 



+ 



(e(fc, A) ■ kc) {e{ka, A') ■k)j] e'^^^o+k). 




,^ , k^ - k ■ kc ^ 
^ k 



{e{k,X)-e{kG,X')) 



(19) 



Notice that the last term in (|T3D does not contribute to PkqW because of the transversality 
condition of the polarization vectors. The integration in d^k is dominated by the upper limit, 
i.e. k ^ kc and accordingly we can ignore the effect of the terms proportional to kc- In 
addition, for those modes k which are outside the Hubble radius at the end of inflation, we 
have kr] <^ 1. With these simplifications we obtain: 

dkdQ . , _ , , , , , 2 



N, 



kc 



E 

A,A' 



J dr] (V2k ((z^ ' + kT) {e{k, A) ■ e(A^G, A')) 



(20) 



(27r)3/2 2kGV 

Performing the integration in the angular variables and using the definition of the power spec 



trum in (18), we obtain: 



ka 



4(27r 



1 3/2 



3fc, 



G 



.Ml (A- 



n-1 



j dr] (z^ ' + kj^) 



(21) 



Finally, we will estimate the time integral. The behaviour of scales that reenter the horizon 
during the radiation dominated era is oscillatory with a decaying amplitude |15|, therefore, 
there is no long-time contribution to the integral that could spoil the perturbative method. 
Thus, for simplicity we will assume that the function JF vanishes for rj > 1/k, and accordingly 
we estimate, |/ drj [ij^ ' + kj^)\^ ~ Our final expression for the occupation number is: 



4(27r)3/M| 
Skcik 



c) 



kc 



dkk 



n— 1 



4(27r)3/2A^ 



^max 



3 n 



kc k^-' 



(22) 



3 Magnetic field generation 



The energy density stored in a magnetic field mode Bk with wavenumber k is given by: 

dpB 



Pb{uj) = u- 



duo 



(23) 
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Figure 1.- \og{B'^'^^/l G) as a function of log{k,nax / kc) ■ The continuous line corresponds to the 
scale-invariant Harrison-Zeldovich spectrum with n ~ 1, the dashed line to n = 0.8 and the dotted 
line to n — 1.25. The dashed horizontal line represents the weakest galactic dynamo seed field 
limit corresponding to a flat universe with cosmological constant and h = 0.5. 



with u = k/a the physical wavenumber. In terms of the occupation number, it reads; 

Pb{uj) = oo^Nk. (24) 

From (^), we can obtain the strength of the field at decoupling on a coherence scale corre- 
sponding to kc ~ 10^3^ GeV as: 

_ V2(a;^ ) iV,^ _ ^^^^^ (25) 



In Fig.l we have plotted the strength of the magnetic field generated as a function of 
the comoving cut-off frequency kmax for different values of the spectral index n. Notice that 
the results are in general too weak to explain the observed fields without any amplification. 
However, for certain values of the cosmological parameters, the produced fields could act as 
seeds for a galactic dynamo. 

We see that the spectrum of magnetic fields produced by this mechanism is thermal Bk ~ 
/c^/^, in the low-momentum region. We can then compare this spectrum with that corresponding 
to the thermal background radiation with a temperature T^ec — 0.26 eV present at decoupling 
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time. The energy density in photons with comoving wavenumer kc at decouphng is given by 
pRi^^a) - k%Tdec/al^^. Thus we find: 

Pr{^g) _ adecTdec _ ^^36 ( J^cY ^26) 



Pb{,i^g) Nkakc ' \k, 



max . 



From this expression we see that the magnetic field energy density will dominate over the 
background thermal radiation whenever log{kmax/kc) > 36/n, i.e. for example, for n = 1 this 
implies kmax > 10^^ GeV. 

The cut-off frequency k^ax cannot be easily determined in general, since it depends on the 
specific mechanism that generates the perturbations and also on the evolution of the universe 
during reheating and thermalization. However we can estimate typical values in some particular 
regime. In the case in which metric perturbations are generated by infiation, it is natural to 
expect, kmax ^ o-iHi, as commented before. Thus, let us take the simplest chaotic infiation 
model with potential = Xcp^/A with A ~ 10"^^ fixed by COBE. In this model the 

Hubble parameter during infiation is Hi ~ 10^^ GeV. Owing to the uncertainties commented 
before, we will let the reheating temperature TptH be a free parameter. After reheating the 
universe evolution is adiabatic aj/adec ~ Tdec/TjiH, and we can calculate the cutoff frequency 
as kmax/kc ~ CLjHj/kc ~ adecTdecHi/{TRHkc), which yields k^ax/kc ~ 10"^^ GeV /Trh- Here 
we have assumed that the inflation-radiation transition takes place in a few inflaton oscillations 



[pO|] (see also ||21|). Comparing with Fig. 1 we see that with these simple estimations for 
the A0^/4 model, the ampliflcation could be above the requirements of the galactic dynamo if 
Trh < 10^ GeV for n = 1.25. 



4 Conductivity effects 

In the previous discussion we have assumed that electric conductivity of the universe played 
no role in the generation of the magnetic flelds. Although neglecting conductivity is a good 
approximation during inflation, it is not during the reheating or radiation periods. As com- 
mented before, the copious production of particles during reheating produces the growing of 
the conductivity which becomes very high during radiation This implies that the mag- 
netic flelds produced in the inflation-radiation transition will evolve conserving magnetic flux 
Pb ~ a^^- However, it has been recently showed that the growth of conductivity during 



reheating could affect the evolution of the EM modes. The effects of conductivity can be taken 
into account in a phenomenological way by introducing a current source Jj = cTcaA- in (|TT]). 



This approach is only valid at sufficiently large scales |]T3|] and in general the rigorous treat- 
ment would requiere to solve the set of coupled EM-matter flelds equations (Vlasov equations) 
which is beyond the scope of this paper. In general, ac{x,ri) is a time and space dependent 
function, and it has been shown that it grows exponentially during parametric resonance in a 
non-equilibrium plasma in QED [^. The calculation of the actual function crc(x, r/) in our case 
would be rather involved and model dependent, since it would require information about the 
reheating and thermalization processes. For that reason we will not take any particular model. 
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but we will do a general discussion of the possible effects in different cases. We will also assume 
for simplicity that all the Fourier modes of the conductivity have the same time evolution, i.e. 
ac{k,r]) = T,{ri)ak during reheating. 



In previous works |ll2| , |T3[, the conductivity was considered as an homogeneous field, adri) 
and because of the form of the current source term, its effect was the damping of each EM 
mode. In our production mechanism, the inhomogeneities are able to mix different EM modes 
and for this reason we need to have information about the complete spectrum and not only 
about the large-scale components. In addition, our analysis is perturbative and therefore we 
can only describe the initial stages in which the conductivity is still small. Let us then consider 
the modified spatial equation to first order in the metric perturbations and the conductivity: 

(2$' - aa,)Af^' + 9,4'^' - 4'^" + 2V<I> ■ V^f ^ - 2V$ ■ dj^^^ + V^A^^ + 4$VMf ^ = (27) 
Following a similar analysis we find the same equation (|T2|), but with the new current: 




2 k 
+ (6(fc, X)-q)^k + q, r^) ^e-^" - ($(^ + ^)e-'='') q, j (2^ 



In ||T^ the following lower limit on the (homogeneous) conductivity is obtained aj ac ~ ai Hj/a, 
with a the fine structure constant, i.e. aj ac ~ kmax/o^ and this implies that the conductivity 

term will dominate in J^''^{kG,Ti) for k <^ kmax- However, as commented before, the dominant 
contribution to the EM amplification comes from the high-frequency modes, i.e. k ~ k^ax- 
In such case, the importance of the conductivity term is determined by the ratio (Xk/Ck when 
k — »■ kmax- Thus, if the conductivity is almost homogeneous, its spectrum will decline at short 
scales and we expect its contribution to the EM field evolution to be negligible. In the opposite 
case, the analysis would become much more involved, since the above magnetohydrodynamical 
approximation would break down and the full set of microscopic equations would be needed. 
In any case, this simple analysis shows that the possible damping effects mainly affect those 
modes with ~ fcc and that, depending on the actual conductivity spectrum, the rest of modes 
could be less severely affected than in other models. 



5 Conclusions 

In this work we have studied the production of photons in the presence of an inhomogeneous 
gravitational background. We have shown how the breaking of conformal invariance induced 
by the evolution of metric perturbations in the inflation-radiation transition is able to produce 
particles, and we have related the occupation number with the scalar metric perturbations 
power spectrum. 

We have considered the possibility that this mechanism could have had some relevance in 
the problem of galactic magnetic fields and we have concluded that the total amplification is 
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several orders of magnitude below the observed strengths. However, for certain values of the 
cosmological parameters and with the assistance of the dynamo mechanism, the amplification 
could be compatible with the current (low-redshift) galactic observations. We have also con- 
sidered the effect of conductivity in a phenomenological way and show that although it could 
affect the evolution of EM modes, in some cases and depending on the particular form of the 
spectrum, the effects could be small. 

The mechanism studied in this work only relies on the existence of a primordial spectrum of 
metric perturbations, described by the scalar spectral index n and a possible cutoff frequency 
kmaxi which are the only free parameters in the model. Therefore, as a by-product we get that 
magnetic fields could also provide useful information about the metric perturbations spectrum 
and in particular about its small-scales region. 

Acknowledgements: This work has been partially supported by the Ministerio de Edu- 
cacion y Ciencia (Spain) (CICYT AEN 97-1693). 
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